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^ . Abstract 

Wideband spectrum sensing is becoming increasingly important to cognitive radio (CR) systems 
for exploiting spectral opportunities. This paper introduces a novel multi-rate sub-Nyquist spectrum 
sensing (MS^) system that implements cooperative wideband spectrum sensing in a CR network. MS'^ 



can detect the wideband spectrum using partial measurements without reconstructing the full frequency 
spectrum. Sub-Nyquist sampling rates are adopted in sampling channels for wrapping the frequency 
spectrum onto itself. This significantly reduces sensing requirements of CR. The effects of sub-Nyquist 
sampling are considered, and the performance of multi-channel sub-Nyquist samplings is analyzed. To 

(N 

. improve its detection performance, sub-Nyquist sampling rates are chosen to be different such that the 



numbers of samples are consecutive prime numbers. Furthermore, when the received signals at CRs are 
faded or shadowed, the performance of MS^ is analytically evaluated. Numerical results show that the 
proposed system can significantly enhance the wideband spectrum sensing performance while requiring 
low computational and implementation complexities. 
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I. Introduction 

The radio frequency (RF) spectrum is a scarce natural resource, currently regulated by gov- 
ernment agencies. Under the current policy, the primary user (PU) of a particular spectral band 
has exclusive rights to the licensed spectrum. With the proliferation of wireless services, the 
demands for the RF spectrum are continually increasing. On the other hand, it has been reported 
that localized temporal and geographic spectrum utilization efficiency is extremely low. For 
example, it has been reported that the maximal occupancy of the spectrum between 30 MHz 
and 3 GHz is only 13.1% and its average occupancy is 5.2% in New York City [1]. The spectral 
under- utilization can be addressed by allowing secondary users to access a licensed band when 
the PU is absent. Cognitive radio (CR) has become one promising solution for realizing this 

goal m, 0. 

A crucial requirement of CRs is that they must rapidly fill spectral holes without causing 
harmful interference to the PUs. This ability is dependent upon spectrum sensing, which is 
considered as one of the most critical components in a CR system. In a multipath or shadow 
fading environment, the signal-to-noise ratio (SNR) of the primary signal as received at CRs 
can be severely degraded, which will not only lead to unreliable spectrum sensing results, but 
will also reduce the capacity of the CR network due to the decreased data transmission time per 
frame. In such a scenario, cooperative spectrum sensing could increase the reliability of spectrum 
sensing by exploiting spatial diversity. In our previous work flU, [|51, centralized cooperative 
spectrum sensing frameworks have been developed for improving the reliability of spectrum 
sensing. However, these studies only considered narrowband spectrum sensing techniques, the 
extension to wideband cooperative spectrum sensing requires yet a different approach. 

To exploit more spectral opportunities over a large range of frequencies (e.g., 10 kHz ~ 10 
GHz), a CR system needs some essential components, i.e., wideband antenna, wideband RF 
front end, and high speed analog-to-digital converter (ADC). Yoon et al. (SI have shown that 
the —10 dB bandwidth of the newly designed antenna can be 14.2 GHz. Hao and Hong fTl 
designed a compact highly selective wideband bandpass filter with a bandwidth of 13.2 GHz. 
In ||8]], Bevilacqua and Niknejad designed a wideband CMOS low-noise amplifier with 
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the bandwidth of approximately 10 GHz. In contrast, the development of ADC technology 
is relatively behind. To the best of our knowledge, when we require an ADC to have a high 
resolution and a reasonable power consumption, the achievable sampling rate of the current ADC 
is 3.6 Gsps dUl. Obviously, ADC becomes a bottleneck in such a wideband system. Even if there 
exists ADC with more than 20 Gsps sampling rate, the real-time digital signal processing of 20 
Gb/s of data could be very expensive. 

In previous work, Quan et al. [[TOll . [[TTIl proposed a multiband joint detection (MJD) approach 
that can sense the primary signal over a wide frequency range. It has been shown that MJD 
has superior performance for multiband spectrum sensing. In [fT2| , Tian and Giannakis studied 
a wavelet detection approach, which could adapt parameters to a dynamic wideband spectrum. 
Furthermore, they cleverly introduced compressed sensing (CS) theory to implement wideband 
spectrum sensing by using sub-Nyquist sampling techniques in the classic paper [fT3l . Later on, 
the CS-based approach has attracted many talented-researchers' attention in [[T4l - [|24l owning 
to its advantage of using fewer samples closer to the information rate, rather than the inverse 
of the bandwidth, to perform wideband spectrum sensing. In [14], Tian et al. studied cyclic 
spectrum sensing techniques with high robustness against sampling rate reduction and noise 
uncertainty. In [fTSl , Zeng et al. proposed a distributed CS-based spectrum sensing approach 
for cooperative multihop CR networks. In our previous work [|25l . [|26ll . to save system energy, 
adaptive CS-based spectrum sensing approaches were proposed that could find the best spectral 
recovery with high confidence. Unfortunately, using CS-based approaches, the spectral recovery 
may cause high computational complexity, leading to a high spectrum sensing overhead due to 
the restricted computational resources in CRs. 

In this paper, we introduce a multi-rate sub-Nyquist spectrum sensing (MS'^) approach for 
cooperative wideband spectrum sensing in a CR network. Because the spectral occupancy is 
low, sub-Nyquist sampling is induced in each sampling channel to wrap the sparse spectrum 
occupancy map onto itself. The sensing requirements are therefore significantly reduced. We 
then analyze the effects caused by sub-Nyquist sampling, and represent the test statistic using a 
reduced data set obtained from multi-channel sub-Nyquist sampling. Furthermore, we propose to 
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use different sampling rates in different sampling channels for improving the spectrum sensing 
performance. Specifically, in the same observation time, the number of samples in multiple sam- 
pling channels are chosen as different consecutive prime numbers. In addition, the performance 
of MS'^ for combining faded or shadowed signals is analyzed, and the closed-form bounds for 
the average probabilities of false alarm and detection are derived. The key advantage of MS^ is 
that the wideband spectrum can be detected directly from a few sub-Nyquist samples without 
spectral recovery. Compared to the existing spectrum sensing methods, MS^ can achieve better 
wideband spectrum sensing performance with a relatively lower implementation complexity. 

The rest of the paper is organized as follows. Section |ll] introduces the signal model. In Section 
Uni we propose the wideband spectrum sensing approach, i.e., MS"^. The performance analysis of 
MS^ for combining faded signals is given in Section |IVl Section |V] presents simulation results, 
and conclusions are given in Section |Vll 



Consider that all CRs keep quiet during the spectrum sensing interval as enforced by protocols, 
e.g., at the medium access control (MAC) layer [fTOl . Therefore, the observed spectral energy 
arises only from PUs and background noise. The bandwidth of the signal as received at CRs 
is W (Hertz). Over an observation time T, if the sampling rate f {f > 2W) is adopted to 
sample the received signal, a sequence of Nyquist samples will be obtained with the length of 
JN = fT. This sequence is then divided into J equal-length segments where denotes the 
number of Nyquist samples per segment (both J and are chosen to be natural numbers). If 
we use Xc,i{t) (t E [0, T]) to represent the continuous-time signal received at CR i, after Nyquist 
sampling, the sampled signal can be denoted by Xi[n\ = Xc,i{n/ f), n = 0, 1, ■ ■ ■ , JN — 1. At 
CR i, the sampled signal of segment j (j e [I, J]) can be written as 



11. Preliminary 




Xe,(n//), n = (j - 1)N, (j - l)Ar + 1, . . . , jAT - 1 



Otherwise. 



(1) 
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The discrete Fourier transform (DFT) spectrum of the sampled signal of segment j is given by 

N-l 

XiAk] = Yl a:,,Ne--'2"'"/^, = 0, 1, ■ ■ ■ , AT - 1 (2) 

n=0 

where j = \l —\. We model spectrum sensing on a frequency bin A; as a binary hypothesis test, 
i.e., ?/o,fc (absence of PUs) and Hi, a,, (presence of PUs) [fTTI : 

V r/i J ^O''^ 

Xij[A;] = <^ (3) 

where Zij[k] is complex additive white Gaussian noise (AWGN) with zero mean and variance 
Sfi^, i.e., Zij[k] ~ CJ\f{0, Sf 1^), Hij[k] denotes the discrete frequency response between the PU 
and CR i, Sij[k] is assumed to be a deterministic signal sent by the PU on the frequency bin k, 
and Hi denotes the spectral support such that = {A;|PU presents at Xj For simplicity, 
in the rest of the paper, we assume that the noise variance of the DFT spectrum is normalized 
to be 1 . The observation time T is chosen to be smaller than the channel coherence time so that 
the magnitude of Hij[k] remains constant within T for one CR, i.e., constant regarding 
the segment number j. 

Because an energy detector does not require any prior information about the transmitted 
primary signal while having lower complexity than other spectrum sensing approaches ^Tf], 
we consider the energy detection approach in this paper. The received signal energy can be 

calculated as j 

Ei[k] = J2\X^Ak]f, k = 0,l,--- ,N~1. (4) 

The decision rule for energy detection approach is then given by 

Ei[k] I Afc, A; = 0,l,---,iV-l (5) 

where is the detection threshold for the frequency bin k. Here, it is noteworthy to empha- 
size that, after Fourier transform, the energy detection is done on each frequency bin in the 
frequency domain. Thus, the noise in high frequencies should not affect the energy detection in 
low frequencies, and vice versa. The benefit of frequency-domain energy detection is that the 
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detection performance depends on the SNR on a single frequency bin, regardless of the noise in 
the other frequencies (e.g., high frequency noise due to wideband sensing). To be specific, the 
signal energy on frequency bin k can be modeled by [|27ll 



k 



EM ^{ ' (6) 

where '-fi[k] = '^(l-^'W'^'WI ) demotes the SNR on the frequency bin k at CR i, X2J denotes central 
chi-square distribution, and X2j(27i[^]) denotes non-central chi-square distribution. Both of these 
distributions have 2 J degrees of freedom and 2'yi\k] denotes a non-centrality parameter. Here, 
the noise has variance 5f i., measured bandwidth and noise temperature r„ = where Kb 
denotes the Boltzmann constant. The probabilities of false alarm and detection are given by fl27l 

ViJ ^) 

Pi,,k = FT{E,[k] > Xk\no,k) = 4x7) 

Pd,i,k = 'Pr{Ei[k] > \k\ni,k) = Qj [V^lWl V^) (8) 

where r(a) denotes the gamma function, T(a, x) denotes the upper incomplete gamma function, 

and Quia, x) is the generalized Marcum Q-function defined by Quia, x) = -^rr t^e ~ Iu-iiat)dt 

in which /t,(a) is the f-th order modified Bessel function of the first kind. 

in. MULTI-RATE SUB-NYQUIST SPECTRUM SENSING 

It is difficult to realize wideband spectrum sensing, because it requires a high speed ADC for 
Nyquist rate sampling. We will now present an MS^ system using multiple low-rate samplers to 
implement wideband spectrum sensing in a CR network. 

A. System Description 

Consider that there are v synchronized CRs collaborating for wideband spectrum sensing, and 
the fusion center (FC) is one of the CRs which has either greater computational resources or 
longer battery life than other CRs. Due to low spectral occupancy [13], the received signals at 
CRs are often sparse in the frequency domain. Here, we assume that the Nyquist DFT spectrum, 
Xi j G C^, is s-sparse (s ^ N), which means that only the largest s out of N components cannot 
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be ignored. The spectral sparsity level, i.e., s, can be obtained from either sparsity estimation [[T6l 
or system initialization (e.g., by long term spectral measurements). As shown in Fig. \T\ MS^ 
consists of several CRs, each of which has one wideband filter, one low-rate sampler, and a fast 
Fourier transform (FFT) device. The wideband filters are set to have bandwidth of W. MS^ can 
be described as follows: 

1) The FC allocates different sub-Nyquist sampling rates to different CRs. 

2) CRs perform sub-Nyquist samplings in the observation time T. 

3) The sub-Nyquist DFT spectrum is calculated by using sub-Nyquist samples and FFT 
device^ 

4) The signal energy vectors are formed by using the sub-Nyquist DFT spectrum. 

5) The CRs transmit these signal energy vectors to the FC by using a dedicated common 
control channel in a band licensed to the CR network 

6) The received data from all CRs is fused in the FC to form a test statistic. 

7) The FC chooses the detection threshold and performs binary hypothesis tests. 

8) The FC shares the detection results with all CRs. 

B. Sub-Nyquist Sampling and Data Combining 

At CR i, we use sub-Nyquist rate fi {fi < 2W < /) to sample the continuous-time signal 
Xc,ii^)- The sampled signal can be denoted by yi[n] = Xc,i{n/ fi), n = 0,1, ■ ■ ■ , J Mi — 1 where 
J Mi = fiT and Mj is assumed to be a natural number. The sampled signal is then divided into 
J equal-length segments. The segment j (j G [I, J]) can be written as 

, ^ { x,^i{n/fi), n={j- 1)M,, (j - l)Mi + 1, . . . , jM, - 1 

= \ (9) 
I 0, Otherwise. 

The DFT spectrum of the sampled signal of segment j (j E [I, J]) can be given by 

YiA^] = Yl ?/^jMe"'''^'""/*'% m = 0, 1, ■ ■ ■ , M, - 1 (10) 

n=0 

'jointly considering wideband spectram sensing and spectrum reuse in CRs, we use FFT devices for distinguisliing different 
frequencies in order to reuse some un-occupied frequencies. Here, tlie use of FFT will cause additional complexity of 
0(M log AI) and memory storage increment, if M denotes the number of FFT points. 
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With the aid of Poisson summation formula [|29l . the DFT spectrum of sub-Nyquist samples can 
be represented by the DFT spectrum of Nyquist samples (as proved in Appendix lAl): 



N 



^X,j[m + /Mi], m = 0,1,- ■■ ,Mi - 1. 



(11) 



/=— oo 



According to (E]) and (fTTT ). the spectral support of the sub-Nyquist spectrum Yij can be given 
by 



{m\m= |A;|^„d{M,)' ^ ^ ^^i- 



(12) 



One risk caused by sub-Nyquist sampling is the signal overlap in ^^^^[m]. However, when we 
choose parameters in JN = fT such that N ^ s and let the sub-Nyquist sampling rate satisfy 
Mi ~ 0{\/N), the probability of signal overlap is very small (as proved in Appendix |B]). In 
such a scenario, we concentrate on considering two cases: no signal on m and one signal on m. 
In the latter case, only a single / is active in (fTTl) . and the other terms in the summation of (fTTT) 
can be modeled as noise by using ([3]). Thus, the following equation holds from (|3]) and (fTTl) : 



J H = ^X^j [m + I Mi] + ^Y1 J ["^ + ^^^] ' rn + ^Mieni 



(13) 



where / is an unknown integer within [0,N/Mi — 1]. Furthermore, using (|3]) and (fTJt . we can 
model the DFT spectrum of sub-Nyquist samples by 

CAfl0H,,[k]S,,[k],^Lk), ken. 
where 5^ ^ ^ is the noise variance of sub-Nyquist DFT spectrum, and can be given by using (fTTl) 




mod (Mi) 



(14) 



ir2 



' N 
M 




5 1 e 



i.k 



6. 



i,k 



(15) 



N \_M^ 

No. of sums Scaling of Yi^j 

where [-^] (the smallest integer not less than -p-) denotes the number of summations in (fTTT) . 
The signal energy of sub-Nyquist DFT spectrum in each CR node is then calculated by 



, m = 0,l,--- ,M,-1. 



(16) 
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X2j5 k ^ Qi 

xij{2f%[k]), ken,. 



(17) 



We note that, due to the sub-Nyquist sampUng, the noise will be folded from the whole 
bandwidth onto all signals of interest as shown in (fl3]) . As a result, comparing (Il7l) with 
©, we find that the received SNR in the sub-Nyquist sampling channel i will degrade from 
7j to ^7i. This SNR degradation depends on the ratio between the number of samples at 
the sub-Nyquist rate and the number of samples at the Nyquist rate (i.e., ^). 

In MS"^, the signal energy vectors at CRs will then be collected at the FC. Finally, we form 
a test statistic by 



In Fig. |2l we give an illustration of the above test statistic for a practical ASTC DTV signal. 
To test whether the PU is present or not, we adopt the following decision rule: 



Let n^^i denote a set of aliased frequencies (i.e., false frequencies appear as mirror images of 
the original frequencies around the sub-Nyquist sampling frequency), and i7u,j represent a set 
of unaffected/unoccupied frequencies: 




(18) 




(19) 




(20) 



(21) 



Using (fTTI) . we can model the test statistic of (fTSl) as 



ken 



■u 




(22) 
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where ^lu = n^^^ilu.i, = U^=iJ^A,i, ^ = ^i=i^i, T = {Afn = \k\ mod iM,),'m G k ^ 
denotes the set of CRs who have aliased frequency on the frequency bin k, and |T| = p denotes 
the cardinality of the set T (equivalently the number of CRs that have aliased frequencies on 
the frequency bin k). 

In (|22|) . k E and k E ^1^ are two extreme cases under the hypothesis 'Ho,k- The former 
case denotes there is no aliased frequency on the frequency bin k, while the latter case represents 
there are maximum number of aliased frequencies (i.e., p) on the frequency bin k. Thus, the 
former one is the best case while the latter one is the worst case for signal detection under the 
hypothesis ?^o,fc- The probability of false alarm on the frequency bin k can therefore be bounded 
by using (|22l) 



< Ptk < Q 



2 '^'"^ 

-J2Ma^[k],^/X~k\. (23) 

iex 



T{Jv) 

We note that the problem of minimizing the probability of false alarm can be transformed to 
minimize the parameter p, which depends on several factors, e.g., the sampling rates of CRs. 
Using the same sub-Nyquist sampling rates in MS'^ is not recommended as it could lead to 
p = V, resulting in the maximum of the probability of false alarm. As we will see in the 
following subsection, the parameter p can be minimized by using different sampling rates at 
CRs. 

C. Multi-rate Sub-Nyquist Spectrum Sensing 

To improve the detection performance of sub-Nyquist sampling system in the preceding 
subsection, we should analyze the influence of sampling rates. Firstly, we consider the case 
of spectral sparsity level s = 1, which means that only one frequency bin ki E is occupied 
by the PU. 

Lemma 1: If the numbers of samples in multiple CRs, i.e.. Mi, M2, My, are different 
primes, and meet the requirement of 

M,Mj>N, \fi^jE[l,v] (24) 
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then two or more CRs cannot have mirrored frequencies in the same frequency bin. 
The proof of Lemma 1 is given in Appendix O 

Secondly, considering the spectral sparsity level s > 2, we find that, if the conditions in 
Lemma 1 are satisfied, the parameter p in (l22l) is bounded by s. It is because only one CR can 
map the original frequency bin kj E f],; to the aliased frequency in ^Ia, and the cardinality of the 
spectral support is s. Therefore, we obtain the detection performance of MS'^ as Theorem [T] 

Theorem 1: In MS^, if the numbers of samples in multiple CRs, i.e., Mi,M2, ■ ■ ■ ,My, are 
different consecutive primes, and meet the requirement of MiMj > N, ^ i ^ j E [l,v], using 
the decision rule of (fT9l ) the probabilities of false alarm and detection have the following bounds: 



Proof: Using (l23l) and the bound |T| = p < s, (l25l) follows. Furthermore, when the energy 
of one spectral component in maps to another spectral component in the probability of 
detection will increase. Thus, the inequality of (|26l) holds. □ 

Remark 1: It can be seen from Theorem [H that the sampling rates in MS^ can be much lower 
than the Nyquist rate because of Mj ~ 0{\/N). By (l25l) we note that the probability of false 
alarm increases when the spectral sparsity s increases. In addition, the higher average sampling 
rate will lead to better detection performance. This is because the probability of signal overlap in 
the aliased spectrum can be reduced with a larger Mj in each sampling channel as our discussions 
in Section IIII-B[ By (f25h and (l26l) . we can see that using more sampling channels (i.e., v), the 
detection performance can be improved. It should be emphasized that there is no closed-form 
expression for the probabilities in Theorem [T] This is because the number of CRs that have 
aliased frequencies on the frequency bin k cannot be predicted. Moreover, we note that the 
upper and lower bounds in Theorem [1] can be easily computed because the Marcum-Q function 
can be efficiently computed using power series expansions [|30l . Under the Neyman-Pearson 
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12 

criterion, we should design a test with the constraint of Pf ^ < a. In such a scenario, we must 
let the upper bound of (|25T ) to be a and solve the detection threshold from the inverse of the 
Marcum-Q function. It has been shown in [31 J that the detection threshold can be calculated 
with low computational complexity. In addition, to calculate the detection threshold, the noise 
power is required to be known at the FC. 

IV. Combination of Faded Signals 

As shown in Section [nil the combining procedure in the FC is to sum up all non-faded signals 
at CRs and make final decisions. In this section, we investigate the combination of faded signals 
at CRs using the same approach as shown in (fTSl ). We assume that the received primary signals 
at different CRs are independent and identically distributed (i.i.d.), and are faded subject to either 
Rayleigh or log-normal distribution. 

Solving the distribution of the sum of weighted independent random variables in (fTST ) is not 
trivial. Hence, we use the sum of uniformly weighted random variables to approximate the sum 
of different weighted random variables in Theorem [1] 

i=l iex 

where M is the average Mj over multiple CRs, ip = 7v = XliLi 7*' 7s = XlleT^"* 7j- 
note that the above approximation accuracy mainly depends on , where smaller 

corresponds to more accurate approximation. Since Mi, M2, are chosen to be v different 
consecutive prime numbers and the distance between primes could be very small compared to 
A^, the parameter will approach to zero as increases. Thus, the above approximation 

has little impact on the final result. 
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A. Rayleigh Distribution 

If the magnitudes of received signals at different CRs follow Rayleigh distribution, then the 
SNRs will follow exponential distribution. Hence, 7v and 7s follow Gamma distributions: 

/(7v) = =i^e-T, 7,>0, /(7s) = =i^e-^, 7s > (28) 

where 7 = E(i^fL) denotes average SNR over multiple CRs, and /(■) denotes a generic 
probability density function (PDF) of its argument. 

The average probabilities of false alarm and detection for MS^ are often solved by averaging 
Pi,k in (|25l) and Pd,fc in (l26l) over all possible SNRs, respectively. 

Theorem 2: If the magnitudes of received signals at different CRs follow Rayleigh distribution. 



the average probabilities of false alarm (Pf,fc) and detection (Pd,fc) in MS'^ will have the following 
bounds 



r{jv) 



where 0(x, Jv, tp, 7, A) is defined as 



<Pf,fc<e(s,Jt;,^,7W,Afc) (29) 



Pi,k>&{v,Jv,^,^[k],Xk) (30) 



e-fi + ^V^'fc" /^y£(!i±M) (31) 

^-1^+9] 2^^n+.-iy^-^^J Tin + Jv) ^^^^ 



n=0 

bl 



in which C„ denotes the binomial coefficient, i.e., C„ = m. 

1 ' ' a a\(b—a)] 

The proof of Theorem |2] is given in Appendix |D] 

Remark 2: From Theorem |2l we can see that < < 1, because the term G [0, 1] and 
the remaining terms can be simplified to 1 . In addition, it can be proved that is a monotonically 
increasing function with respect to ip,j, and x. Therefore, both probabilities will either increase 
or remain the same when the average sampling rate and the average SNR increase, more sampling 
channels will lead to a higher probability of detection, and the average probability of false alarm 
can be reduced with smaller s. 

Remark 3: Because (|3T| ) contains infinite sums, its computational complexity is directly related 
to the number of computed terms that are required in order to obtain a specific accuracy. As the 
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number of computed terms, i.e., P, varies, the truncation error can be written as 



n=P 



where the inequality of (l33l) holds because < 1, and (|34l) is obtained by using the binomial 

expansion. It can be shown that (|3T1) converges very quickly. For example, in order to achieve 
double-precision accuracy, only P = 30 ~ 40 calculated terms are required; therefore the bounds 
are tractable. To solve for the detection threshold A^, we could use the lower bound on Pf ^ in 
(|30] ). This is because the lower bound can approximate Pf ^ very well as analyzed in Appendix ID] 
and also verified by Fig. |3] 



B. Log-normal Distribution 

The strength of the transmitted primary signal is also affected by shadowing from buildings, 
hills, and other objects. A common model is that the received power fluctuates with a log-normal 
distribution. In such a scenario, the PDF of the SNR at CR i, i.e., /(7i), is given by 

e /_(101og.o(7.)-7,."' 

where ^ = 10/ln(10), and cTj (dB) denotes the standard deviation of 101og^Q7j at CR i. Note 
that the PDF in (|35l) can be closely approximated by a Wald distribution [|27l . 



/(7.) = ^^exp|- ^-"'^t7 I. 7.>0 (35) 



where 6i = E(7j) denotes the expectation of 7i, and rji is the shape parameter for CR i. Via the 
method of moments, the parameters r^j, 6i and 7^, a, are related as follows: 

9, = exp h-± + ^\, n, = ^ . (37) 



{ 2ey' " exp(f)-l 
In the proposed system, the condition ^ = ^^^^ = b (constant) can be satisfied. Thus, 7s 
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and 7v will also follow the Wald distribution [|33l . The PDFs of 7s and 7v are given by 



where r] and 9 denote the averages of r/j and 9i, respectively. 

Theorem 3: If the magnitudes of received signals at different CRs follow log-normal distribu- 
tion, the average probabilities of false alarm (Pt;^) and detection (Pd,fc) in MS'^ will be bounded 
as 

iV> A(t;,Jt;,^,Afc,e[A;],r/[A;]) (41) 
where A(x, Jf , V^, A, 6, rf) is defined by 



1 



in which K„_i (a) denotes the modified Bessel function of the second kind with order n — ^ 

"2 ^ 

The proof of Theorem |3] is given in Appendix |El 

Remark 4: Because (|42l) contains infinite sums, the truncation error T\(P) must be considered. 
Similar to (|33] ). the truncation error can be written as 



TT n! " '2 

n=P 

1 



vr ^ n! "-^1^ ) 

It can be shown that (l43l) decreases to zero very quickly as P increases, therefore the bounds 
in Theorem |3] are easy to compute. 

V. Simulation results 

In our simulations, we assume that the CRs are organized as shown in Fig. [T] and adopt the 
following configurations unless otherwise stated. We use the wideband analog signal model in 
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Il34l and thus the received signal Xc,i(t) at CR i has the form: 

XoAt) = Yl \Hi,i\VEiBi ■ sine {Bi{t - A)) • cos {27r fi{t - A)) + z{t) (44) 
1=1 

where sinc(x) = A denotes a random time offset, z{t) is AWGN, i.e., z{t) ~ A/'(0, 1), 

El is the transmit power at PU, and Hn denotes the discrete frequency response between the 
PU and CR i in subband /. We generate v = 22 independent channels according to the fading 
environment, and regenerate them for next observation time. The received signal Xc^t) consists 
of A'^^ = 6 non-overlapping subbands. The l-th subband is in the frequency range of [/; — -y, 
fi + ^], where the bandwidth 5; = 1 ~ 10 MHz and fi denotes the center frequency. The center 
frequency of the subband / is randomly located within [-y, W — ^] (i.e., /; G [y, W — ^]), where 
the overall signal bandwidth W = 10 GHz. If the wideband signal were sampled at the Nyquist 
rate / = 2W for T = 20 /iS, after segment division with J = 5, the number of Nyquist samples 
per segment would be = 80, 000; thus, using FFT-based approach, the frequency resolution 
is Y/j = 0-25 MHz. In MS^, the received signal is sampled by using different sub-Nyquist 
rates at different CRs. To be specific, the numbers of samples in multiple CRs are chosen by 
using Theorem 1 and we choose the first prime Mi ^ ay/N (a > 1) and its i; — 1 neighboring 
and consecutive primes. The spectral observations are obtained by applying an FFT to these 
sub-Nyquist samples in each channel. Then the signal energy is calculated in the spectral 
domain using (Il6l) . and the energy vectors are transmitted from the CRs to the FC using 
dedicated common control channels. These channels are assumed to be Rayleigh block 
fading channels (constant channel gains over one time block) corrupted by circularly 
symmetric complex Gaussian noise with zero mean and unit variance. In addition, we 
consider that the channel power gain of the common control channel is normalized to unit 
and the average SNR as received at the FC is 15 dB. In the FC, we form the test statistic by 
using (fTSi) . We define the compression rate as the ratio between the number of samples at the 
sub-Nyquist rate and the number of samples at the Nyquist rate, i.e., ^ where M denotes the 
average number of sub-Nyquist samples at CRs. Spectrum sensing results are obtained by using 
the decision rule (fT9l ) and varying the detection threshold A^. 
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In Fig. [31 we verify the theoretical results in (|25])-(|26]), (|29l)-(l30l), and (|40l)-(|4l]) by comparing 
them with the simulated results. It shows that the lower bound on the probability of false alarm 
can tightly predict the simulated results while the upper bound seems relatively loose. It is 
because that the assumption (i.e., all s components in the Nyquist DFT spectrum will be mapped 
to the same location when the signal is sub-Nyquist sampled) for deriving the upper bound can 
rarely occur. Fig. |3] also illustrates that the lower bound on the probability of detection can 
successfully predict the trend of simulated results. Comparing the faded signal cases with the 
non-faded signal case, it is found that, when combining faded signals, the probability of detection 
declines more slow than that combining non-faded signals. This is more obvious for the case of 
combining signals following log-normal distribution as shown in Fig. [3tc). 

Fig. HI shows the receiver operating characteristic (ROC) curves of MS^ when combining non- 
faded and faded signals. When the average SNR as received at CRs is 5 dB, the performance of 
MS^ combining faded signals is roughly the same as that of combining non-faded signals. This 
is because the strength of the signal is mostly masked by the noise. In contrast, the detection 
performance of MS'^ combining non-faded signals outperforms that of combining faded signals 
when SNR=10 dB. In addition, it is seen that the performance of MS^ combining log-normal 
shadowed signals is the poorest. Nonetheless, even for log-normal shadowed signals, MS^ has a 
probability of nearly 90% for detecting the presence of PUs when the probability of false alarm 
is 10%, with the compression rate of ^ = 0.0219. To investigate the influence of s and SNR, 
we use Fig. \5\ to show the performance of MS'^ when the received signals are faded according 
to Rayleigh distribution with different values of s (proportional to the number of subbands). We 
see that, as the number of subbands decreases, the detection performance improves for the same 
SNR. The performance improvement of MS^ stems from that, for a fixed number of sampling 
channels, decreasing s makes it easier to distinguish the occupied frequencies from the aliased 
frequencies as discussions in Section IIII-C[ 

Fig. [6ta) depicts the influence of the standard deviation when the MS^ system combines 
log-normal shadowed signals. It can be seen that a larger standard deviation will lead to worse 
detection performance for the MS'^ system. It is because a larger a is equivalent to a longer tail in 



February 7, 2013 



DRAFT 



18 

the log-normal distribution, thus making the detection more difficult. In Fig. ^h), we compare 
the performance of MS^ with that of Nyquist systems. In the Nyquist system type I, each 
CR is given an orthogonal subband (wideband spectrum is divided into several equal-length 
subbands) to sense using Nyquist rate, while their decisions are sent back to the FC. In the 
Nyquist system type II, we assume that each CR must sense all wideband spectrum non- 
cooperatively, thus requiring multiple standard ADCs in each node to cover all wideband 
spectrum. After signal sampling, all measurements are sent back to the FC, where equal 
gain combining approach is adopted to fuse data and then energy detection is used for 
spectrum sensing. Fig. |6](b) shows that the proposed system has superior performance to 
the Nyquist system type I, but inferior performance to the Nyquist system type II. The poor 
performance of the Nyquist system type I mainly results from the lack of spatial diversity 
gain. In the Nyquist system type I, each subband is only sensed by one CR as each CR is 
given an orthogonal subband to sense, which cannot take advantage of spatial diversity. 
In contrast, both the proposed system and the Nyquist system type II are monitoring each 
subband using several CRs, thus taking advantage of spatial diversity. It can also be seen 
that the Nyquist system type II has marginal performance gain over the proposed system, 
however, at the expense of much higher implementation complexity as discussed below. 

In Table HI we compare the implementation complexity of MS'^ with that of the Nyquist 
systems, when the received signals at different CRs are faded according to Rayleigh distribution. 
Here, we consider the comparison metric: the number of same- sampling-rate ADCs for achieving 
Pd > 90% and Pf < 10%, because practical CRs often have requirements on the probabilities 
of detection and false alarm to secure the performance of both CRs and PUs. We can see that, 
when there exist 10 CRs, MS'^ requires each CR equipped with a single ADC with an average 
sampling rate of 957.54 MHz; thus, the whole CR network only requires 10 low-rate ADCs. In 
contrast, the Nyquist system type I requires 21 ADCs in total, because of 21 x 957.54 MHz^ 20 
GHz for covering 10 GHz spectrum based on Nyquist sampling theorem. In the Nyquist system 
type II, 210 ADCs (with the average sampling rate 957.54 MHz) will be required because each 
CR will require 21 ADCs similar to the Nyquist system type I. Thus, the system complexity 
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of MS^ is approximately half of that of the Nyquist system type I and much less than that 
of the Nyquist system type II. 

In Fig. |71 we choose the CS-based system in lfT3l as a benchmark system due to its high 
impact and outstanding performance. The comparison between the proposed MS^ system and 
the benchmark system is provided. We assume that t> = 22 CRs are collaborating for wideband 
spectrum sensing in both systems, in order to increase the reliability of spectrum sensing by 
exploiting spatial diversity. We can see from Fig. |7Ia) that MS'^ outperforms the CS-based system 
for every compression rate. In Fig.|7tb), it is seen that, compared with the benchmark system, 
MS^ has better compression capability. Using MS^, the probability of successful sensing 
becomes larger than 90% when the compression rate ^ > 0.023. In contrast, the benchmark 
system can achieve the probability of successful sensing 90% only when the compression 
rate ^ > 0.045. Furthermore, as shown in Table III we can find that the computational 
complexity of MS^ is 0{N log N) due to the energy detection with FFT operations, rather than 
O {N{M + log A^)) in the CS-based system, where M is usually much larger than log A^. The 
complexity of the CS-based system is caused by both the matrix multiplication operations and the 
FFT operations for spectral recovery. To sum up, with the same computational resources, MS'^ 
has a relatively smaller spectrum sensing overhead than the CS-based system, not only because 
of the better compression capability (less data transmission results in shorter transmission time), 
but also due to the lower computational complexity. 

VI. Conclusions 

In this paper, we have presented a novel system, i.e., MS'^, for wideband spectrum sensing 
in CR networks. MS^ can relax the wideband spectrum sensing requirements of CRs due to 
its capability of sub-Nyquist sampling. It has been shown that, using sub-Nyquist samples, the 
wideband spectrum can be sensed in a collaborative manner without spectral recovery, leading 
to a high energy-efficiency and a low spectrum sensing overhead. Moreover, we have derived 
closed-form bounds for the performance of MS'^ when combining faded or shadowed signals. 

Simulation results have verified the derived bounds on the probabilities of false alarm and 
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detection. It has also been shown that using partial measurements, MS'^ has superior performance 
even under low SNR scenarios. The performance of MS^ improves as either the number of CRs 
or the average sampling rate increases. Compared to the existing wideband spectrum sensing 
methods, MS'^ not only provides computation and memory savings, but also reduces the hardware 
acquisition requirements and the energy costs at CRs. 

Appendix A 

Relationship between Nyquist DFT Spectrum and Sub-Nyquist DFT Spectrum 
Using the Poisson summation formula [|29l , dD, and (fTOl ). we obtain: 

/, J2 ^c,(^ + = Yl = Yl = Y,{w) (45) 

zez neZ n=0 

where Xc,i{w) = Xc,i(t)e~-'^™*(it. Similar to (HSj, by using and (0), we can obtain: 

N-l 

/5^Xe,(«; + //) = J^x.Me-^^™'^ = 5^x,[n]e-^2™" = X,{w). (46) 

«ez neZ n=0 

As the received signal is bandlimited and / > 2W, Xi(w) = fXc^i{w) holds for w G [— 
Substituting it to (1451) . we obtain Yi{w) = j X]/^-oo + ^ discrete form, we end up 

with: ^ 

M 

Yi[m\ = J2 ^^i'^ + ^^il^ m = 0,l,--- ,Mi-l. (47) 
Appendix B 

Probability of Signal Overlap at Sub-Nyquist Sampling 

As s spectral components are distributed over the frequency bins of 0, 1, ■ ■ ■ , — 1, the 
probability of the frequency bin k belonging to the spectral support Q is P = Pi(k E ^l) = 
Let q denote the number of spectral components overlapped on the frequency bin m, using (fTT)) 
the probability of no signal overlap is given by 

Pr(g < 2) = Pr(g = 0) + Pr(g = 1) = (l-P)^^^ + (^^'^\p{l-P)^^.^-^ (48) 
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where \^~\ denotes the number of summations in (fTTT ). Substituting P = into (|48] ) while 



choosing sub-Nyquist sampling rate in MS such that Mi = y N, we obtain 



P,(, < 2) ^ [^t^4r f jXzf)^ (49, 



y Mi \ N J N-s 
It can be tested that Pr(g < 2) approaches to 1 when choosing such that N ^ s. Thus the 
probability of signal overlap approaches to zero under the condition we choose. 

Appendix C 
Proof of Lemma 1 

Let Mi and M, denote the number of samples at CRs i and j, respectively. Using (fT2] ) and 
(|20l) . we can represent the aliased frequencies projected from ki E Q by 

Qi = \ki\^od{M,) + lMi = ki-hMi + lMi, h^l (50) 
9j = \ki\^„diM,) + lMj = ki-hMj+lMj, h^l (51) 

where integers h and h are quotients from modulo operations, and l — hE [— \^~\ +1, — 1], 
I — h e [— \jT-] + 1, \^] — 1], in which \^] gives the smallest integer not less than 

Avoiding gi = gj is equivalent to avoiding (/ — h)Mi = (/ — h)Mj. If Mj and M, are 
different primes, the condition max(|/ — h\) < Mj (i.e., — 1 < Mj) will satisfy this. After 
simplification, the condition MiMj > N is obtained. Moreover, if it holds for any two CRs, the 
case for more than two CRs will also hold. 

Appendix D 
Proof of Theorem [2] 

If the received signals at CRs are Rayleigh faded, the lower bound on the average probability 
of false alarm will remain as it is independent of the SNR. Using (|25] ). (|271) , and (|28l) , the upper 
bound on the average probability of false alarm can be calculated by 

P^"' = r Qjv (^/ins, =|^e- T rf^,. (52) 

Jo ^ ' 7 r(s) 
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Rewriting the Marcum Q-function by using (4.74) in [|35l and (8.352-2) in [[36ll . we obtain: 



ni L (n+ Jv) 

Substituting (|53l) into (1521) . we can rewrite (|52] ) as 

7 ^ '^'(s - l)!r(n + J?;) 
Calculating the integral by using (3.351-3) in [|36l . we end up with 

Similarly, we can obtain the lower bound on the average probability of detection. 

Appendix E 
Proof of Theorem [3] 

If the received signals are shadowed according to log-normal distribution, the lower bound 
on Pf in (|40l) will remain. By (l38l) . the upper bound on the probability of false alarm can be 
given by 

- v^) y|or-e.p (-^) ... (56, 

u 

Substituting (1531 ) into (1561) , we calculate Ptyc as 

^^^^ =v^g ^^.,rUj.) V o ^''^ 

Using (3.471-9) in ||36ll for calculating the integral in (|57T ), we obtain 

^^•'^ "VV^ n!r(n + J.;) [\] l^JO^ ) I^V J' ^^^^ 

Likewise, the lower bound on the average probability of detection can be approximated. 
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Fig. 1. Block diagram of multi-rate sub-Nyquist spectrum sensing (MS'^) system. 
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Fig. 2. Illustration of the proposed system for sensing practical ASTC DTV signal: (a) Spectrum of an ASTC DTV 
signal in Nyquist sampling system, and (b) Test statistic in the proposed system. Here, we use ASTC DTV signal 
WAS_3_27_06022000_REF and assume the number of CRs v = 22. The number of Nyquist samples is AT = 80, 000, and 
the numbers of sub-Nyquist samples in MS'' system are consecutive primes Mi — 1613, M2 = 1619, • • ■ , AI22 ~ 1783. The 
average compression rate is calculated to ^ = 2.12%. 
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Fig. 3. Comparisons of simulation results and analytical results for the probabilities of false alarm and detection when MS"' 
combining (a)non-faded signals, (b) Rayleigh faded signals, and (c) Log-normal shadowed signals with the received SNR= 5 
dB (at CRs) and o- = 4 dB. 
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Fig. 4. Receiver operating characteristic curves of MS^ for combining non-faded signals or faded signals when the compression 
rate = 0.0219 and the number of segments J = 5. The wideband signal is observed by 22 CRs at different sampling rates 
(the average sampling rate is 448.68 MHz). 
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Fig. 5. The performance of MS^ for combining Rayleigh faded signals with v — 22 and ^ = 0.0219, when the received 
SNR at CRs and the number of subbands change. 
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Fig. 6. Receiver operating characteristic curves for combining log-normal faded signals: (a) MS"^ when the standard deviation, 
i.e., a, and the average SNR at CRs vary, (b) comparison between MS^ and Nyquist systems when ct = 5 dB and the average 
SNR as received at CRs is 5 dB. In simulations, the wideband signal is sampled at different sampling rates by 22 ADCs with 
the average sampling rate of 448.68 MHz, and the compression rate is ^ = 0.0219 and J — 5. 



TABLE I 

Implementation complexity comparison of MS^ and the Nyquist systems when the received signals are 

FADED according TO RAYLEIGH DISTRIBUTION WITH TEN DECIBEL RECEIVED SNR AT CRS. 



Number of CRs for 

Wideband Spectrum Sensing (v) 


10 


20 


30 


40 


Required Number of ADCs 
in the Proposed MS^ System 


10 


20 


30 


40 


Required Number of ADCs in 
the Nyquist System Type I 


21 


40 


58 


74 


Required Number of ADCs in 
the Nyquist System Type II 


21x10 


40x20 


58x30 


74x40 


Average Sampling Rate of ADCs 
in the Above Systems (/ in MHz) 


957.54 


513.08 


350.34 


276.77 


Sample Reduction Rate I 

^ Total samples in MS"' 

^ Total samples in Nyquist type I ^ 


47.62% 


50% 


51.72% 


54.05% 


Sample Reduction Rate II 

^ Total samples in MS"^ 

^ Total samples in Nyquist type II ^ 


4.76% 


2.5% 


1.72% 


1.35% 
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Fig. 7. Comparison between MS'^ and CS-based system [13|: (a) the probability of detection when the probability of false 
alarm is set to 10%, and (b) the probability of successful sensing which is defined as the probability of achieving both Pj > 90% 
and Pf < 10%. In simulations, the average SNR as received at CRs is 10 dB and the number of CRs is u = 22. 



TABLE II 

Comparisons of wideband spectrum sensing techniques. 



Approach 


Compression 


ADC/DSP 


Implementation 


Computational 


Capability 


Type 


Complexity 


Complexity 


Wavelet detection 


X 


Nyquist 


low 


C(A^logA^) 


Multiband joint detection 


X 


Nyquist 


high 


C(A^logA^) 


CS-based detection 


V 


sub-Nyquist 


medium 


0{N{M + log N)) 


Proposed system 


V 


sub-Nyquist 


low 


0{NhgN) 
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